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I- IN IBOPaC TION 



Cne of the U.S. Coast Guard's missions is maintaining 
aids to navigation in the waters of the United States. 
Ihese aids reguire periodic servicing to ensure they are on 
station and they are showing their proper watch characteris- 
tics. The responsibility of maintaining many of these aids 
falls on the Coast Guard's fleet of buoy tenders. A buoy- 

tender in its area of responsibility may have as many as 200 
aids which it is required to maintain. 

The buoytender has scheduled ATON runs (Aids TO 
Navigation) and emergency outages it must handle during a 
fiscal year. An emergency outage occurs when an aid is 
reported by a mariner to be showing improper characteris- 
tics or is missing from its station. These outages are 
handled shortly after they are reported. Several times 
during the year the buoytender plans an ATON run to handle 
the periodic servicing of aids, to ensure that they maintain 
their proper watch characteristics. 

Prior to a scheduled run the tender will receive SANDS 
forms from the unit's district office for each of the aids 
which are due for service or relief. SANDS is a database 
where information is kept on a district's navigational aids. 
The district's aids to navigation branch identifies aids 
which are due for service or relief and notifies the tender. 
The district provides a list of aids which need to be 
serviced during the scheduled run and of work needed on each 
aid. The tender then makes plans accordingly to schedule 
and complete the work. 

The number of aids serviced on a particular scheduled 
ATON run may be as few as 5 (a day of local work) to as many 
as 70 (an extended four-week trip servicing the west coast 
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of Alaska from the Aleutian Islands to Point Hope) . Ihe 
maintenance of the floating aids and fixed structures 
entails the travel of the tender from aid to aid. The tend- 
er's Operations Officer usually plans the initial route to 
he travelled to service the aids due for routine maintenance 
and presents his proposal to the Commanding Officer, who may 
then modify this proposal. 

There are several factors which go into the decision 
process to produce the desired route. Two of these factors 
are (1) which aids are scheduled for service, and (2) the 
distances between the aids. The desired route is usually 
the route of shortest distance which visits all the sched- 
uled aids and returns the tender to her point of origin. 
This problem is a classic problem which operations analysts 
normally call 'the traveling salesman problem.' Present 
route selection methods are (1) traditional or previously 
followed routes and (2) routes created by Operations 
Department personnel sitting down with chart and rule and 
selecting a route. 

It is the intent of this thesis to propose and implement 
a computer assisted approach as an alternative method for 
the buoytender route selection problem. The succeeding 
chapters of this paper contain a discussion of the trav- 
eling salesman problem and how it relates to the buoytender 
problem, a brief survey of possible solution methods to the 
problem, a description of the solution process of the 
selected method, and the results of the route selection 
program with conclusions. 
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II- HATORE OF THE PR^LEM 



In this chapter the principles of the traveling salesman 
problem will be presented and it will be shown how the buoy- 
tender problem can be described as a traveling salesman 
problem. The discussion will include a general description, 
definitions about solution methods, and the solution 
reguirements of the problem. Further, a description of the 
computing hardware on which the selected solution method is 
programmed is provided. This provides an idea of the capa- 
bilities of the hardware since the hardware has a bearing on 
the selection of the solution method. 

A. DESCRIPTION OF THE TRAVELING SALESMAN PROBLEB 

A brief description of the traveling salesman problem 
(or TSP) is as follows. A salesman has n cities he must 
visit. The salesman starts in one of the cities and must 
travel so that he enters and departs each of the other 
cities only once. Upon visiting the last city the salesman 
will return to the city from which he originally started his 
travels. The desired solution to the problem is the route 
the salesman can traverse which is of the shortest total 
distance possible. 

It is usually assumed that a TSP can be associated with 
some n by n distance matrix. The elements of the matrix, 
d(' , are distances from city i to city j, (i = 1,...,n), 
(j= 1 , . . , ,n) ,and where djXsoo, i=1,...,n. The traveling 

salesman is not allowed to leave city i and to return to 
city i in his tour. 

A TSP can be associated with a graph, also. A city 
which the salesman must visit is represented by a node in 
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this graph. There is an arc (i,j) in the graph with length 

for each node. 

A distance matrix may be classified as sparse or dense. 
A sparse matrix will have a majority of its entries equal to 
infinity. The corresponding graph is sparse since it has 
few edges. A dense matrix is a matrix where a majority if 
not all its entries, with the exception of the dit ’s, are 
less than infinity. The number of possible permutations of 
cities for a dense matrix is of order (n-1)!. Computation 
time and memory required for solution by total enumeration 
grows rapidly with the size of n due to the large number of 
possible solutions, each permutation of cities being a 
feasible tour and a possible solution. 

The distance matrix used by the buoytender for route 
selection will be a dense matrix. With the large number of 
permutations or possible solutions to the problem, some 
means other than total enumeration must be used to solve the 
buoytender problem. 

The traveling salesman problem may be broken into two 
classes, symmetric and asymmetric. In a symmetric TSP the 

while in an asymmetric TSP this need 
not be true. Further, a TSP may or may not be required to 
satisfy the triangle inequality. Consider three cities i, 
j, k. The triangle inequality states that 

dij < dik + dkj for all i/j^k^i. 

The conditions which make a TSP asymmetric, the unequal 
distances between cities, may cause the problem to violate 
the triangle inequality. 

The TSP, both symmetric and asymmetric, belongs to the 
N? Complete class of problems [Ref. 1,2]. This has one 
major implication with respect to the complexity of solving 
the TSP; no polynomial-time algorithms are known or seem 
likely tc be devised for exact solution of the TSP. Only 
exponential-time exact algorithms are known. 



There have been many algorithms proposed foe the solu- 
tion of the TSP, each algorithm showing some advantages and 
disadvantages to its method of attacking and solving the 
problem. Several of the published algorithms are tailored 
to a particular size of TSP. The TSP may be categorized by 
size as either small or large- A reference could not be 
found in the literature which would give a hard numerical 
value by which the size of a TSP could be judged, but a 
small TSP is generally considered to be less then 15 cities. 
A large ISP is often considered to be in the neighborhood of 
45 or more cities. There are solution methods which solve 
the TSP exactly for small up to large problems. There are 
some methods which can solve some large problems exactly, 
but as a rule large problems are solved with heuristics or 
approximate methods. 

The circumstances of the buoytender problem satisfy the 
conditions of a class of the traveling salesman problem. 
The buoytender problem can be characterized as a symmetric 
TSP, i.e., the distance traveled between any two aids is the 
same regardless cf direction travelled by the tender. The 
size of buoytender problem will vary with the number of aids 
scheduled to be serviced on a specific ATON run; therefore 
the solution method for the buoytender problem will have to 
handle both small and large problems. 

The triangle inequality as applied to the buoytender 
problem holds true, tut how it holds needs some explanation. 
With all distances positive, the triangle inequality states 
that given three points A, B, and C, the distance from A to 
C must be less than or equal to the distance from A to B 
plus the distance from B to C (see Figure 2.1) . In the 
buoytender 's problem the straight-line distance from A to C 
may be less than (A to B) + (B to C) , but some situations 
will require that A to C equal (A to B) + (B to C) (see 
Figure 2.1), 
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Figure 2.1 Triangle Ineguality / Obstructed A to C Path 




The buoytender may be unable to take the direct route from A 
to C due to shoal water or an obstruction. To get from A to 
C the tender will have to take the routes A to B and B to C. 
This is not an unusual situation for buoytender s since the 
aids being serviced often mark shoal water or obstructions. 



B. DESCHIPTIOH OF TEE SOLOTICB REQOIREHENTS 



There are many proposed solution methods to the TSP and 
they may be divided into two groups; exact methods and 
approximate or heuristic methods. An exact method or algo- 
rithm has the property that upon termination of the algo- 
rithm, the user will have the best possible or global 
optimua solution. Heuristics or approximate methods, on the 
other hand, terminate with a local optimum solution or the 
best solution found thus far. This local optimum solution 
may be the global optimum but usually there is no proof that 
the global optimum has been found. Heuristics or approxi- 
mate methods solve by checking a subset of the numerous 
possible solutions and terminate when the subset has been 
searched for the best answer it contains, a specified time 
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The C3 Standard Terminal is a standalone 16 bit micro- 
computer. The terminal comes in two basic configurations, 
the Integrated Hork Station (IWS) and the Application Work 
Station (AWS) . The IWS is a master or standalone unit 
equipped with a hard disk and 8 inch floppy mass storage 
system, a 300/1200 baud modem, a printer, and one of three 
RAM memory configurations - 256k, 384k, or 512k. The IWS 

may be programmmed in BASIC, FORTRAN, COBOL, Pascal, and 
Assembly. The ARS is more limited in its capabilities than 
the IWS. The AWS is limited to 256k RAM and must be 
networked to a IWS station to access the peripherals. With 
a program in RAM the AWS will operate as a standalone 
computer, but does not have the power or computational speed 
of the IWS. 

The solution method selected for programming will have 
to perform within the limitations of the C3’s memory. The 
C3*s RAM memory will be a critical factor in selection of 
the solution method. Since the computer is a purchased 
standalone and not a leased system, CPU time should not be a 
critical factor. It is still desired, though, to keep CPU 
time from becoming excessive since the computer is also 
needed for ether work. The computer could be allowed to run 
overnight (off duty hours) for large problems, and therefore 
a solution method requiring more than 12 hours to provide an 
acceptable solution would be excessive. A solution method 
which provides an answer in under an hour could be run 
during duty hours with very little impact on other users, 
particularly if it was run on a slave terminal. 

In the next chapter we take a brief look at some 
suggested solution methods for the TSP and their applica- 
bility to the buoytender problem. 
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III. A BRiEF S0R7EY OF SOLDI ION METHODS FOR FHE ISP 



In this chapter a brief survey of the general solution 
methods for the TSP and some of their associated algorithms 
will be presented. The number of algorithms proposed for 
solving the TSP is extensive. To present all of the avail- 
able algorithms or to discuss them in detail is worthy of a 
dissertation in itself. This survey of solution methods is 
made because a combination of a couple of the methods will 
be used to solve the buoytender problem. This survey will 
be brief in nature presenting the general solution methods 
and briefly describing some of the more popular of published 
algorithms which fall under these solution method headings. i 

The solution methods may be classified under four 
general headings. These headings are Tour Building, Subtour 
Elimination, The One-Tree Formulation, and Tour Improvement. 
These headings describe how the traveling salesman problem 
is approached. 

A. TCDR BUILDING 

Tcur Building solution methods to construct a tour using 
the n nodes and the available arcs in the problem. Tour 
Building contains both exact and approximate algorithms. 
Some of the approaches taken by Tour Building Algorithms to 
solve the TSP are dynamic programming, ’branch and bound' 
methods, and tour construction heuristics. 



^Three excellent articles [Ref. 1,3,4] have been 
published which describe the TSP, discuss general solution 
methods, and review many of the algorithms which sclve it. 
These three articles review the research conducted on the 
TSP from the mid 1950's up to 1983. 
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Held and Karp [Eef. 5] present a dynamic programming 
approach which is suitable for small problems and can be 
modified to an approximate method for larger proDlems. Ihe 
dynamic programming approach rapidly consumes available 
memory as the size of the problem increases. Due to memory 
limitations, the algorithm Held and Karp presented is exact 
only up to 13 cities. They also present a dynamic program- 
ming method which approximates the solution for larger prob- 
lems. The approximation is done by partitioning the large 
problem into several smaller problems or subproblems, each 
subproblem containing 13 or fewer nodes. Each partitioned 
subproblem is solved for its optimal tour resulting in a set 
of subtours for the original problem. The subtours are then 
broken and linked to each other to provide a tour for the 
original large problem. 

Little et a 1. [Bef. 6] present a branch and bound 
method. The algorithm branches on whether a particular arc 
is included or excluded from the tour being constructed. 
This branching partitions the solution set containing all 
tours into subsets. Each arc in the total graph can be 
associated with some subset of tours, and the idea is to 
find the subset containing the optimal tour. For each 
subset of tours, a lower bound is computed for the tours 
within the subset. The tour is constructed as arcs with 
favorable Icwer bounds are included in the tour. The algo- 
rithm continues branching and computing lower bounds until a 
subset is found which contains one tour and this tour's 
distance is less than or equal to the lower bounds of the 
other subsets of tours. This tour will be the optimal solu- 
tion to the TSP. The memory required for this technique may 
be extensive and computing time grows exponentially with the 
size cf the problem. 

Many tour construction heuristics have been proposed. 
Eosenkrantz, Stearns, and Lewis [Eef. 2] and Golden et al. 
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[Bef. 4] review several of these heuristics. The two most 
popular forms of these heuristics are the nearest neighbor 
algorithm and the insertion technigues. 

The nearest neighbor algorithm operates exactly as its 
name implies, A starting node is chosen and its nearest 
neighboring node is selected as the next node in the tour 
being constructed. The algorithm then iteratively selects a 
node not yet in the tour which is nearest to the previously 
selected node to enter the tour. This selection continues 
until all nodes have been selected and then the last node 
selected is connected to the first node to complete the 
tour. 

The insertion technigues include the nearest, farthest, 
cheapest, and arbitrary insertion algorithms. Insertion 
technigues begin with a starting node and select the nearest 
node to create a subtour. The algorithm iterates selecting 
a node which is nearest to or farthest from any other node 
in the subtour, or arbitrarily selects the next node to 
enter the tour. This selected node is then inserted into 
the subtcur wherever it provides the least increase (or in 
the case of the farthest, the greatest increase) in distance 
in the new subtour. These heuristics generally provide a 
suboptimal solution to the TSP, but they have the advantage 
of being very simple to implement and have very fast solu- 
tion times with minimal memory consumption, 

B. SOBTOOE ELIHINATICH 

The Subtour Elimination solution method [Ref, 1] attacks 
the TSP by solving an n by n assignment problem with the 
added constraints that the final solution must contain a 
cycle and the cycle cannot be of size n-1 or less. The 
number of integer constraints necessary to eliminate all 
subtours or cycles of size n-1 or less is enormous. 
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Therefore, the initial solution is generally found using a 
’relaxed* form of the assignment problem. The relaxed 
assignment problem emits the subtour constraints in its 
formulation. 

Mathematically stated the relaxed assignment problem is 



Min 
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The solution to the relaxed assignment problem provides 
the initial lower bound on the optimal tour length. If the 
optimal assignment solution is a tour then it is an optimal 
solution to the TSP. If the assignment solution is not a 
tour then there exist subtours which must be removed until a 
single tour exists. If a subtour exists with k arcs in the 
subtour, then there are k possible subproblems to which the 
problem may branch. These k subproblems each have an addi- 

tional constraint, the constraint for subproblem i 
(i=1,...,k) being the exclusion of arc i from the problem to 
eliminate the subtour. Instead of branching into k subprob- 
lems the approach is to branch into two subproblems. An arc 
is selected from the subtour and the exclusion of this arc 
becomes a new constraint in one of the subproblems while 
inclusion of the arc in the tour becomes a new constraint 
for the alternate subproblem. The modified assignment 
problem is solved again for each of these new constraints 
and provides the lower bound for its respective branch. In 
this way only the necessary subtour elimination constraints 
are added as needed to the assignment problem, rather than 
attempting to add all possible subtour constraints to the 
initial assignment problem. Ideally the branch solution 
subset with the lowest bound is then selected to be tested 
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to see if it is a tour, and, if it is not, the prcblem 
tranches again on the lowest bound. This branching 
continues until all subtours have been eliminated and a tour 
is found which has a distance less than or equal to the 
lower bounds of the ether solution subsets. 

There are two basic methods for searching the tree 
created by the branch and bound process. One method is a 
breadth-first search. As the program is branching and 
computing bounds, a tree is created with leaves or terminal 
nodes, A breadth- first search requires storage for the 
solutions of each leaf and a search through these leaves for 
the next subproblem on which to branch. This requires 
extensive, usually exponential, storage. Branching at the 
lowest bound requires a breadth-first search. The second 
method for searching the tree is a depth-first search. The 
depth-first search is a more localized approach to handling 
the search. The depth-first search need only store the 
immediate solution while searching a branch. If during the 
search it is found that no further branching is possible 
from the current node, the process must return to a previous 
leaf to continue the search for the optimal tour. The 
depth-first search backtracks and recreates the previous 
solution at that leaf. The depth-first search does not 
require exponential storage like the breadth-first search, 
but it usually requires more computation time. 

The subtour elimination method is described as being 
exact. The subtour elimination method usually requires less 
computation time than the branch and bound approach of 
little et al. since the bounds obtained with the assignment 
relaxation are tighter. 
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C. TKE CNE-TREE FORBOLATION 



The one-tree formulation was introduced and developed by 
Held and Karp [Ref. 7,8]. The general method for the solu- 
tion of the TSP begins with the construction of a minimum- 
weight ’cne-tree.* If the one-tree is not a tour then an 
integer linear program may be used to obtain the final solu- 
tion. 

To assist in the description of this solution method 
several terms need defining. First a tree is defined as a 
connected graph without cycles. A minimum- weight spanning 
tree is a tree with ainimum total weight on its edges. A 
one- tree is a tree to which one edge has been added yielding 
exactly one cycle. A minimum-weight one- tree is a one-tree 
with minimum total weight on its edges. Finding a minimum- 
weight one-tree is a relaxed version of a TSP since the 
soluticn to a TSP is a min imum- weight one-tree having every 
vertex of degree two. 

In the solution method proposed by Held and Karp 
[Ref. 7] the algorithm begins by constructing a minimum- 
weight one-tree. A minimum-weight one-tree can be can be 
found by first constructing a minimum- weight spanning tree 
through nodes 2 to n and then adding to the graph the two 
arcs of least weight from node 1 [Ref. 9], A minimum- 
weight one-tree may also be constructed by constructing a 
minimum- weight spanning tree through nodes 1 to n and then 
add to the graph the next minimum-weight edge not yet used 
from the distance matrix. If the one-tree constructed is a 
tour (i.e. the one-tree’s vertices are all of degree two) 
then the TSP is solved. Otherwise, the one- tree must be 
converted into a tour. 

Integer linear programming is used to obtain an optimal 
tour if the one-tree formulation does not satisfy a tour. 
Held and Karp introduce the concept of a ’gap* function and 
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then use a special integer linear program to minimize this 
function. in essence there is a transformation on the one- 
tree variables which allows the integer linear program to 
search for the optimal tour. Held and Karp suggested two 
methods for optimizing the integer linear program, an ascent 
method and a branch and bound method with the ascent method 
embedded in it. Other improvements to the Held and Karp 
method were suggested by Hansen, Krarup [Ref. 10] and Houck, 
Picard, Vemuganti [Ref. 11], Bazarra and Goode [Bef. 12] 
did further work on Held and Karp's approach using optimiza- 
tion of a lagrangian dual in lieu of solving integer linear 
program. They proposed a branch and bound scheme with subg- 
radient optimization of the dual to transform the one-tree 
into the optimal tour. The one-tree method as described is 
exact and may require extensive computation time. 

D. TCOR IMPROVEMENT 

The Tour Improvement method assumes that an arbitrary 
tour is available. The method operates by perturbing or 
exchanging arcs in the tour until a better tour is found. 
The Tour Improvement method terminates when a better tour 
cannot be found. In its simplest form, it is possible for 
this method to check (n-1) 1/2 arc exchanges for the symme- 
tric TSP before terminating (i.e., it investigates all 
available answers) . 

Dantzig, Fulkerson, and Johnson [Ref. 13,14] proposed a 
algorithm which starts with an arbitrary tour and uses 
integer linear programming to improve the tour. The TSP is 
transformed into an integer linear program and solved using 
the simplex method. Dantzig, Fulkerson, and Johnson 
[Ref. 15] mention that the number of constraints needed to 
characterize the problem is ’astronomical. ' Instead of using 
all the constraints, they begin by using a relaxed version 



22 



of the problem and add constraints as needed to the integer 
linear program to maintain feasibility. In their method 
infeasibility appears as ’loop* or subtour. To continue 
solving the problem a constraint is added which removes the 
subtour yet does not eliminate any of the available tours. 
This algorithm requires the addition of constraints during 
the solution process to maintain feasibility. Several 
researchers have suggested improvements on Dantzig, 
Fulkerson, and Johnson’s method for solution of the TSP. 
[Ref. 16,17,18,19] 

Croes [Ref. 20] proposed a heuristic whereby a subset of 
all possible transformations is tested, a transformation 
being the transformation of one tour into another tour. He 
called these transformations ’inversions’ because they 
inverted the sequence of nodes in part of the tour to create 
a possible improvement. Lin [Ref. 21] later describes Croes 
’inversion free tours’ as 2-optimal tours and goes on to 
describe k-optimality where k is some fixed number less than 
n. The general idea is to exchange k arcs iteratively in 
the tour while testing for improvements. Lin and Kernighan 
[Ref. 22] presents a modified k-opt algorithm where k is not 
fixed. In this algorithm k may vary from 2 to n. The 

floating k-opt algorithm is also described by Christofides 
and Eilon. [Ref. 23] 

The advantage of the ’inversions’ or k-opt solution 
method is that for a given problem, the memory needed to 
solve the problem is fixed. There are no constraints which 
must be added to the problem, and the decision rule used is 
simple. The k-opt method is exact, but may be terminated 
before optimality is reached to provide a satisfactory 
(suboptimal) solution. This generally reduces run time and 
still provides a good solution to the problem. The k-opt 
method may also be programmed to require very little user 
interaction; in the best case the user need only enter the 
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distance matrix for tbe nodes desired in the tour. Ihe 
k-opt algorithm may require extensive computation time, 
particularly if during solution k is found to he greater 
than 2 on many of the iterations. 

This chapter has presented a general survey of some the 
solution methods for the TSP and some algorithms associated 
with these solution methods. The next chapter will discuss 
which of these solution methods were employed to solve the 
bouytender problem. The buoytender routing algorithm will 
then be presented and discussed in detail. 
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17. BOOY TENDER ALGORITHM 

This chapter will discuss the algorithm with which we 
propose to solve the buoytender problem. The reasons for 
the algorithm's selection and some background on its devel- 
opment will be presented, together with a detailed descrip- 
tion of how the algorithm addresses the buoytender problem. 



A. ALGORITHM SELECTION 

The algorithm is a combination of two of the previously 
mentioned solution methods in that it combines a Tour 
Building heuristic with a Tour to Tour Improvement 

heuristic. The concept of combining a Tour Building 
heuristic with a 2-cpt or 3-opt heuristic was proposed by 
Golden et al. [Ref. 24] as a relatively fast computational 
solution method which would provide an optimal or near 
optimal tour. Their 'composite' algorithm is the foundation 
on which the buoytender algorithm was developed. 

The heuristics used in the buoytender algorithm are the 
nearest neighbor algorithm and the k-opt algorithm. The 
nearest neighbor algorithm is used to construct an initial 
tour while a version of the k-opt proposed by Lin and 
Kernighan [Ref. 22] is used to improve this initial tour. 
The nearest neighbor and k-opt heuristics were selected for 
use in the buoytender algorithm because very little interac- 
tive work is required, the interactive work is simple, no 
constraints need to be added to maintain feasibility, compu- 
tational time is reasonable, and the algorithm operates with 
a fixed amount of memory. 

The buoytender algorithm as presented requires very 
little user interaction. The nearest neighbor algorithm 
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needs only the entry of the distance matrix for the problem 
and the selection of a starting node. The k-opt algorithm 
needs the same distance matrix and an initial tour. The 
heuristics* decision processes will provide a feasible and 
possibly optimal solution in a reasonable period of time. 
The nearest neighbor algorithm was selected to construct the 
initial tour because it is computationally faster than 
insertion techniques, and provides a good initial solution 
to the problem. [Ref. 25] 

A critical factor in solving the ISP on a microcomputer 
is the memory required by the algorithm used. The amount of 
memory required for solution of a TSP of size n is fixed for 
the nearest neighbor and k-opt heuristics. With known 
memory requirements for the heuristics, a program can be 
developed to fit the Coast Guard Standard Terminal. The 
nearest neighbor algorithm requires an n by n distance 
matrix and an n arrary in which to store the tour as it is 
constructed. The k-cpt algorithm requires the same storage 
as the nearest neighbor plus a n by n decision matrix, and a 
2n array for recording nodes selected for inversion in the 
tour. 

E. PfiOGBAHHlNG THE AIGOEITBH 

The huoytender algorithm was developed as a series of 
modules which were then programmed as subroutines in the 
operating program. The algorithm was programmed in FORTRAN 
for operation on the Coast Guard Standard Terminal. FORTRAN 
was selected because it is efficient and a majority of 
programmers are familiar with it. 

There are some differences between the heuristics in the 
literature and the buoytender algorithm. The nearest 
neighbor algorithm as used by the buoytender algorithm is 
the same as outlined by Rosenkrantz, Stearns, and lewis 



26 



[Eef. 2], The k-opt heuristic used in the buoytender algo- 
rithm is a modified version of the heuristic presented by 
Lin and Kernighan [8ef. 22]. There are three major differ- 
ences between the k-cpt in the bouytender algorithm and the 
k-opt suggested by Lin and Kernighan. First, the Lin and 
Kernighan k-opt heuristic selects all possible nodes which 
will result in a shelter tour before creating the improved 
tour. The buoytender algorithm improves the tour by 
exchanging nodes whenever a favorable selection is found. 
This improvement-as- ycu-go procedure results in a simplifed 
selection process for the next set of nodes to be tested. 
Second, the Lin and Kernighan k-opt heuristic has an added 
facility for 'limited backtracking' used when a particular 
exchange gives an improvement of zero. The backtracking 
procedure searches until a gain greater than zero can be 
found which improves the tour. In the interest of short- 
ening computational time, the buoytender algorithm does not 
'backtrack' when a gain of zero is found, but instead treats 
zero gain as no improvement and continues with the next 
selection. The third major difference between the Lin and 
Kernighan k-opt heuristic and the buoytender algorithm is 
related to what Lin and Kernighan call 'reduction.' Their 
heuristic, after producing several locally optimal tours, 
checks fer arcs which appear in each of the tours. These 
'good' arcs are not allowed to be broken in further computa- 
tions for other locally optimal tours thereby reducing the 
number of links to be checked for improving the tour. VJhile 
decreasing run time, this procedure requires more memory and 
is therefore omitted from the buoytender algorithm. 

The buoytender algorithm can be divided into three basic 
operations. First, the algorithm selects a random sample of 
k nodes, k<n. This random sample of nodes becomes the set 
of starting nodes for the next major operation in the algo- 
rithm, to create k nearest neighbor tours. The third 
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operation is to use the k nearest neighbor tours as initial 
tours for the k-opt improvement. The best of the k-opt 
improvement tours is presented as the solution to the 
problem. 

A further modification was made to the algorithm based 
on some computational results obtained from use of an early 
version of the algorithm. During programming of the buoy- 
tender algorithm several benchmark statistics were obtained 
from the program. Two of these statistics were the initial 
tour distances produced in the nearest neighbor phase, and 
the distances of the final tour solutions produced by the 
k-opt phase of the algorithm. There seemed to be a rela- 
tionship between the tour produced by the nearest neighbor 
phase and the improved tour produce by the k-opt phase. A 
least-squares linear fit of the initail tour distance vs. 
the final tour distance was done for four data sets. In 
each case the fit pr educed a positive slope. This implies 
for a relatively large nearest neighbor tour, the final tour 
from the k-opt phase will be relatively large. A relatively 
small nearest neighbor tour will produce a relatively small 
tour from the k-opt phase. Using this information the buoy- 
tender algorithm was further modified to reduce computation 
time without compromising too much of its ability to produce 
optimal solutions. for a set of k nodes, k<n, the nearest 
neighbor algorithm provides k initial tours. The set of k 
nodes is selected at random from the population of n nodes. 
Then for a set of L nodes, L<k, the k-opt algorithm is run 
to obtain L final tours. The set of L nodes are drawn from 
the k nodes which produced the L shortest initial tours. 
The shortest final tour from the L final tours produced by 
the k-opt phase will be the algorithm's solution to the 
problem. The final version of the algorithm has a greatly 
reduced operating time over the initial version. 
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C. AIGOBITHM DESCRIPTION 



In this section, a detailed description of the buoy- 
tender algorithm is given. In this description references 
will be given to the program subroutines so that the reader 
may associate the portion of the algorithm under discussion 
with its operational counterpart in the program. 

The buoytender algorithm begins by reading a file from 
mass storage containing the number of aids (nodes) in the 
problem, the distance matrix, and the names of the aids 
(subroutine OBTAIN). A program HATUTL.FOE, which assists 
the user in creation of this data, can be found after the 
listing of the buoytender route selection program, 

ROUTE, FOE. 

After obtaining the data the algorithm selects a random 
sample of nodes, based on the size of n, from which it will 
create tours (subroutine NODE) . The function which selects 
the size of the random sample is 

Sample Size = 5 ♦ f2 ( Log ( n-4 ))\ . 

The function is designed to capture almost the entire popu- 
lation for testing when the problem is small and to capture 
only a small portion of the population being tested when the 
problem is large. This produces a sample size which 
increases the probability of obtaining the global optimum 
when the problem is small yet produces a sample size which 
will maintain reasonable computation times when the problem 
is large. 

After determining the sample size, a random sample of 
nodes is taken. Since there is no random number generator 
for FORTRAN installed on the Coast Suard Standard Terminal, 
a pseudo-random number generator was added to program (func- 
tion RANDOM). The coding for the generator comes from 
Wichmann and Hill [Ref, 26] ; it is described by its 
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authors as an efficient psuedo-random number generator 
having a cycle length of 2.78 x 10^3 . Nodes are sampled at 
random without replacement from the integer population of 1 
to n for use as starting nodes for the nearest neighbor 
phase. 

The algorithm then moves into the nearest neighbor phase 
(subroutine NBR) . A tour is created for each node in the 
random sample, with the node from the sample as the starting 
node of the tour. The algorithm selects the nearest node to 
the starting node for the next entry in the tour construc- 
tion. The algorithm then iterates selecting a node not yet 
in the tour which is nearest to the previously selected node 
as the next entry in the tour. This process continues until 
all nodes have been selected, at which time the last node is 
connected to the starting node to complete the tour, and the 
tour distance is computed. After all the initial tours have 
been constructed they are ordered by their total distances 
from shortest to longest. The five sho rtest tours are 
retained for possible improvement in the k-opt phase. 

A decision matrix is created which designates which 
distance arcs in the distance matrix are presently members 
of an initial tour (subroutine HARKD) . The algorithm then 
proceeds to step through each node of the tour testing it 
for k-optimali ty . 

The tour is 'prepared* at each node to be checked for 
k-optimality (subroutine TRPREP), This process simplifies 
some of the later operations. Here, it is helpful to 
consider the tour as a circle with n postions on its circum- 
ference (see Figure 4.1), and to let positions on the 
circular tour be denoted by p (i) where i = 1,...,n, A node in 
the tour is represented as a position on the circle. TRPREP 
places the starting node in p(1) and places the remaining 
nodes in the tour at p (2) through p (n) where (i) is the 
node’s position in the tour. 
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Figure 4. 1 Circular Representation of a Tour 

Since floating k-cpt algorithm is in essence a series of 
2-opt iterations a description of the simpler 2-opt prin- 
ciple will be presented first. Figure 4.2 shows the 2-opt 
principle in graphic detail. In Figure 4.2 (A), the tour is 
represented in circular form. The tour is now perturbed to 

try and find an improvement. To perturb the tour in the 

2-opt procedure, two arcs will be broken in the existing 
tour and two new arcs will be selected as incoming to form a 
new tcur. This provides a gain function which can be used 
to check for improvement. The gain function is 

G = .E (broken arcs) - 1 (incoming arcs) . 

If G is positive this set of arc exchanges will improve the 

tour - 

Tc perturb the tour, one of the arcs incident to node 
being checked for 2-optimality must be broken or removed 
from the tour. In the buoytender algorithm the node being 
checked for k-optimality will always be in p(1) of the 
circular tour and the arc being broken will always be 
between p(1) and p(2) of the circular tour. An incoming arc 
must then he selected to replace the broKen arc. This 
incoming arc must be selected from those arcs with a 
starting node which is the same node as the end node of the 



3 1 






(E) Feasible Arc Exchange 



'F) Circular Tour 



Figure 4.2 Graphic Description of 2-opt 



previously broken arc. Figure 4.2 (B) shows seme of the 
available selections for the iccoming arc. The incoming arc 
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is in general the shortest available arc, but may be any of 
the available arcs which maximize G. This subtour must be 
broken and another incoming arc must be selected so as to 
reconstitute a tour. Therefore, one of the tour's arcs 
incident to the end node of the first incoming arc must be 
broken. The end node in Figure 4.2 (C) is at p(7). There 
is only one of two possible tour arcs which may be broken to 
allow the incoming arc to become part of the tour. Figure 
4.2 (D) shows the selection of the wrong incident tour arc 
for removal. The breaking of this arc and the subseguent 
reconnection to the node being checked for 2-optimality does 
not result in a tour. Instead, two subtours are created 
giving an infeasible situation. Figure 4.2 (E) shows the 
breaking of the proper incident arc and subseguent reconnec- 
tion tc the original node to create the 'improved tour'. 
Also notice that some of the arc directions in the original 
tour must be reversed to provide a consistent direction of 
movement through the tour. Reversing these arc directions 
is the same as inverting a seguence of nodes in the tour. 
The result is shown in Figure 4.2 (F) as a new 'improved' 
circular tour. 

The 2-opt procedure is carried oat using each node in 
the tour as the starting node checked "or 2-optimality. One 
interation of the 2-opt procedure through all n nodes will 
produce a tour which is 2-optimal. 

The k-opt algorithm is very similar to the 2-opt proce- 
dure described above. Instead of just reconnecting the tour 
as the 2-opt procedure does (p (6) to p(1)), the k-opt checks 
to see if another incoming arc, starting at the end node of 
the last broken arc, can be selected to further improve the 
tour (Figure 4.3). The procedure interates until no further 
improvement can be found. A k-opt, using an example with k 
equal to four, might look like Figure 4.4 The floating k-opt 
algorithm has the ability to vary k during the k-opt 

r 
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Figure 4.3 K-opt Selecting Next Incoming Arc 







solution process. Ihis is more desirable than fixing k in 
the solution process. For example, if k were fixed at four, 
the algorithm can only check four arc exchanges for improve- 
ment but the problem may need a five arc exchange to produce 
the optimal solution. By allowing the k-opt to vary k, the 
probability of the algorithm finding the optimal solution is 
increased. The floating k-opt is bounded between 
2-optimality and n-optimali ty. The k-opt cannot be greater 
than n-optimal because an n-opt exchange creates the optimal 
tour. 
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The k-opt procedure in the buoytender algorithm operates 
by checking each node in a tour for k-optimaiity . The algo- 
rithm, using subroutine TRPEEP, creates the circular tour 
for the node being checked for k-optimali ty. The circular 
tour is prepared so that the initial arc being broken will 
always be between p(1) and p(2). The algorithm checks the 
circular tour at the starting node for k-optimalit y. If no 
improvement can be made breaking the first arc between p (1) 
and p (2) then the circular tour’s direction is reversed 
(subroutine REVTOR) so that the alternate incident arc is 
placed between p(1) and p(2) of the circular tour and this 
tour is checked for k-optim ality. 

The breaking of the arc between p(1) and p(2) requires 
an incoming arc to replace it (subroutine SELCTY) . The 
algorithm selects the five shortest available arcs which 
originate at position two of the circular tour and terminate 
at another position on the circular tour. The arc from p (2) 
to p (3) is not available since it is already in the tour, 
and the arc from p(2) to p (1) is not available since it is 
the arc just broken to improve the tour. For each of the 
selected incoming arcs another arc must be broken to elimi- 
nate a subtour (as in the 2-opt procedure) . It is then 
necessary to reconnect the arcs to reconstitute the tour. 
This results, for each of the selected arcs, a set of two 
incoming arcs and two outgoing arcs. Using the previously 
defined G, each of the five sets of arcs is checked to see 
which maximizes G. If G>0 then the exchange or inversion is 
made (subroutine ADJTCR) , otherwise no exchange is made. 

If an exchange is made the algorithm goes back and 
searches for a new incoming arc from the end point of the 
last broken arc. Subroutine ADJTOR has inverted the tour so 
that the last arc used to reconnect the tour is now between 
p(l) and p (2) of the circular tour. Subroutine SELCTY is 
again used to select the next best available incoming arc as 
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described above. The algorithm interates through selecting 
an incoming arc and inverting the tour as long as improve- 
ment, i.e. positive gain, can be found. 

When improvement can no longer be found the direction of 
the tour is reversed at the last node where improvement was 
found, and the tour is tested again for possible further 
improvement. When nc further improvement can be found and 
the reversed tour has also been tested, the algorithm then 
increments to the next node in the tour and starts the k-opt 
procedure again testing this node for k-opt imality. Each 
node in the tour is used as a starting node in the k-opt 
procedure. When all nodes in the tour have been tested for 
k-optimality the k-opt phase ends and the resulting tour is 
stored. 

After each of the five shortest initial tours have been 
tested for k-optimality, the shortest of the improved k-opt 
tours is selected as the solution to the buoytender problem. 
Figure 4,5 is a flowchart of the buoytender algorithm. 

The next chapter discusses the results of some test 
problems solved by the algorithm and makes concluding 
remarks about the algorithm. 
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Figure 4,5 



Algorith* Flowchart Section 1 
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Figure 4.6 



Algorithm Flowchart Section 2 
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V. PBOGEM RES PIT S AND CONCLUSI ONS 



This chapter will discuss some results obtained from 
using the tuoytender route selection program, including 
results from an actual problem from the U.S. Coast Guard 
Cutter BIACKHAW. All program results were obtained using 
the Coast Guard C3 Standard Terminal. The chapter will 
conclude with some observations regarding the buoytender 
problem, the algorithm, and the test results. 

A. lEOGBAM RESPITS 

To program and test the buoytender algorithm it was 
necessary to obtain problems with known optimal solutions. 
Initially, a very simple and small problem with nine nodes 
was created. This test problem’s optimum tour is the most 
obvious tour obtained by inspection of the problem since the 
nodes are arranged in an almost circular pattern. 
Connecting these nodes following this nearly circular path 
results in the optimum solution. Several alternate tours 
were checked by hand computation and were found to be subop- 
timal. Therefore, the near circular path for this test 
problem is assumed to be optimal. When the program was 
tested with this test problem, the program produced this 
hypothesised optimum solution with a computation time of 13 
seconds. 

With the knowledge that the program appeared to be oper- 
ating correctly on tie initial test problem, further testing 
was necessary to check the program's ability to solve other 
problems. Since the tuoytender route selection problem is a 
symmetric case of the TSP, the literature was searched for 
test problems. Three articles [Ref. 5,20,27] yielded eight 
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symmetric TSP's of various si2es. These eight literature 
problems ranged in size from 5 to 57 nodes with known 
optimal answers. The program, when tested with these prob- 
lems, produced solutions which ranged from the optimal to 
within 3.5 percent of the optimal distance. 

Another problem was created by the random selection of 
70 points on a Euclidean plane. A distance matrix was 
computed for these 70 nodes and this problem was solved by 
the program. The problem’s optimal solution is not known, 
but the problem was run to ensure the program could handle 
the maximum problem for which it was designed, and to 
produce a sample solution time for this maximum sized 
problem. 

With evidence that the program was solving test prob- 
lems, a further test was conducted using an actual buoy- 
tender routing problem, with data provided by the U.S. Coast 
Guard Cutter BLACKHAW (WLB-390), homeported in San 
Fransiscc. The aids to navigation run data was for a trip 
scheduled from 9 April 1984 to 16 April 1984. With the 
experience and historical data this vessel has available 
from servicing these aids for many years, it is hypothesized 
that the route which was scheduled for this ATON run would 
be optimal or nearly optimal. The planned route for the 
ATON run had a distance of 455 nautical miles. When this 
problem was entered and solved by the buoytender program, 
the solution was also 455 nautical miles. Although the 
distances were the same, the scheduled route and the program 
route differed slightly in their sequence for visiting the 
aids. This difference may be attributed to the fact that, 
as stated in Chapter II, the triangle inequality may be 
required to be a strict equality on some of the arcs. This 
difference is possible since the buoytender may, due to the 
strict equality on some arcs, be required to ’pass' a previ- 
ously serviced aid tc complete the route. Figure 5.1 shows 



a five node example of this situation. This suggests there 
may be several sequences of nodes which will provide an 
optimal solution to the buoytender problem. 

TABLE 1 

Tours for BLACKHAB Problem 



Eoute as Scheduled 
Buoy ten d e r 



Route Produced 
by the P rog r am 



1 


7B ISLAND 


1 


YB ISLAND 


2 


HONIARA LWB10A 


20 


MOSS LBB MLA 


3 


AN ISLAND LWE8 


19 


MONTRY BY LB B 


4 


SANTA CROZ ME 


1 8 


MONTRY HBR MB 


5 


P BLANCAS LWB4A 


17 


PT PINOS LHB2 


6 


SAN SIMEON LEB1 


16 


PT CYPRESS LGB6 


7 


VON HELM R LGE4 


5 


P BLANCAS LWB4A 


8 


ESTERO BY GB10E 


6 


SAN SIMEON LBB1 


9 


SOOZA R LGB 


7 


VON HELM R LGB4 


10 


HESTDAHL R LEEl 


15 


ESTERO BY LBB 


1 1 


PT SANLOIS LBE3 


8 


ESTERO GB10E 


12 


LANSING R B 


14 


MORRO BY LBB1 


13 


PT BOCHON LWE2 


13 


PT BUCHON LWB2 


14 


MORRO BY LBB1 


10 


HESTDAHL R LBB1 


15 


ESTERO BY LWE 


9 


SOUZA R LGB 


16 


PT CYPRESS LGE6 


12 


LANSING R B 


17 


PT PINOS LWB2 


1 1 


PT SANLOIS R B 


18 


MON TRY H3R ME 


4 


SANTA CRUZ MB 


19 


MONTRY BY LB E 


3 


AN ISLAND LWB8 


20 


MOSS LBB MLA 


2 


MONTANA LWB10A 


1 


YB ISLAND 


1 


YB ISLAND 





Figure 5. 1 Multiple Optimal Solutions 
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TIBLE 2 

Results of Test Problems 



Problem 

( Re f.) 

Kara (27) 
Tes t (aat 
Bar achet 
3LACKHAW 
Croes (20 
HeldKar 



Karp 

( 27 ? 



Kar g 
Dantzig 
HeldKarp 
Kara (27T 
Max (auth 

♦ 10 05f is 
♦* The di 



& 




S elution 


Tour Dist 


.. (Miles) 


Solution 




Size 


Time (Sec) 


Program 


Optimal 


Quality* 




5 


5 


148 


1 48 


10 0% 


h. I 


9 


13 


22. 95 


22.95 


10 0% 


(27] 


1 1 0 


17 


378 


378 


100% 




20 


50 


455 


455 


♦ ♦ 


) 


20 


52 


254 


246 


10 3. 25% 


(5) 


25 


72 


1711 


1711 


100% 




33 


121 


109 29 


10861 


IOC. 62% 


271 


42 


188 


705 


699 


100. 86% 


(5) 


48 


238 


11847 


11470 


103.28% 




57 


340 


13367 


12956 


103. 17% 


.) 


70 


494 


6791.07 


unk. 


. unk 


optimal; 


105. 28% is 3 


.28% greater than 


optimal. 


stance is 


hypothes ized 


to be optimal. 






Figure 5.2 Least Squares Fit of Size vs. Solution Time 



Table 2 presents a summary of the eleven test problems 
showing problem sizes, solution times, and the quality of 
the solutions produced by the program. All of the run times 
are reasonable by the criteria stated in Chapter II. The 
largest problem, 70 nodes, had a solution time of U94 
seconds, or roughly a little more than eight minutes. 
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Ccmputation times for the eleven problems were plotted 
against problem size, as shown in Figure 5,2 A least-squares 
polynomial fit of the data points is 

Time = ,58h + ,62n + .0914n2 + .00000994n3 
and this curve is also shown in the figure. 

B. CCNCIUSIONS 

The buoytender route selection algorithm provides 
another useful way of scheduling ATOM routes. As a tool for 
the operator to assist in planning the route to be taken, 
the program should produce optimal or near optimal solutions 
for prcblems up to a size of 70 aids. The user may accept 
this route as is or may desire to modify the route based on 
circumst ances requiring human judgement. Appendix A is a 
user’s guide for the matrix utility and buoytender route 
selection programs. 

The buoytender route selection algorithm, with the 
nearest neighbor and k-opt heuristics, provides a quick and 
satisfactory solution to a symmetric traveliag salesman 
problem cn problems up to a size of 70 nodes. The problem 
may have several routes with different sequences for 
visiting the aids, tut the routes may all have the same 
distance. 

It is hoped that this algorithm will be of use to the 
fleet cf U.S. Coast Guard buoytenders. As a decision aid it 
may help the operators of buoytenders obtain optimal or 
nearly optimal routes to service their aids, possibly saving 
time and fuel. 
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APPENDIX A 

USEE’S GUIDE TO PROGBAHS 



The intention of this appendix is to provide instruc- 
tions for the use of the buoytender route selection program 
and utility programs. It should be noted that prior to 
using the buoytender route selection program, the user will 
need to use the matrix utilities program listed at the end 
of this appendix to prepare a data file. 

The buoytender route selection program and the matrix 
utilities program are written in Convergent FORTRAN. 
Convergent FORTRAN is FORTRAN 77 compatible. The programs 
should be typed in as presented and linked to the operating 
system. The Convergent FORTRAN manual has instructions 
concerning compiling and linking. Recommended names for the 
run files are ROUTE. RCN and MATUTL.RUN. Both programs will 
operate cn either a IWS or AWS station with 264 K memory or 
more. The programs run about three time slower on an AWS 
terminal than on an IWS terminal. 

The user will be required to enter text and numerics 
into the programs. All numeric entries are to be integer 
entries (no decimal point) with the exception of the 
distance entries. All distance entries require a decimal 
point entry for proper input. 

MATUTL.RUN is capable of entering problems up to a size 
of 100 aids. ROUTE. RON is designed to handle problems up to 
a size of 70 aids. MATUTL.RUN allows the entry of a large 
matrix of aids and includes a utility which allows the user 
to select a subset of this large matrix for use in the 
program ROUTE. RUN. This allows the user to type in one 
large matrix and then create smaller matrices from the 
larger as needed. This feature will preclude the user from 
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typing in a distance matrix every time he desires to run a 
problem. 

A. OSE CF MATUTl. EOS 

The matrix utilities program is designed for creating 
and editing data files which will be used in the buoytender 
route selection program. A data file will contain the size 
of the problem n, the distance matrix for the problem, and 
the names of the aids to be visited. 

The size of the problem, n, will be the number of aids 
to be visited plus one. The plus one is for the port from 
which the tuoytender will start and end the ATON run. 

To create the distance matrix some extensive plotting is 
required. Most buoytenders have a set of tracklines they 
regularly follow to visit their aids. In most cases these 
existing tracklines will be sufficient to use in the 
distance matrix. To create a distance matrix for entry into 
a data file it is recommended the user use the format shown 
in Figure A.1 . Distances should be computed in nautical 
miles and may range from .0 1 nm to 9999.99 nm. Figure A.1 
is an example of a problem of size nine. 

Since the problem is symmetric, the user need only 
calculate and enter the distances as shown in Figure A. 1 . 
The matrix utilities will automatically complete the rest of 
the matrix entries. The distance you should compute should 
be the shortest possible distance from the row entry to the 
column entry. Using the example in Figure A.1 the entry in 
row 3, column 6 is the shortest possible trackline from aid 
A to aid E. Once the entries for the distance matrix have 
been prepared the user is ready to use MATUTL. RUN. 

MATUTL. EON has some simple menus and some descriptive 
comments to help the user. The menus require a numeric 
response for the selection of an option. The first menu 
presents the utilities available to the user. Figure A.2 . 
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( 1 



1 - create a distance matrix 

2 - extract a small matrix from a larger matrix 

3 - correct a previously stored matrix 

4 - create a hardccpy or a distance matrix 

5 - list names with their row number in the matrix 

6 - exit this pregram 

What is your selection? 



Figure A. 2 MiTOTL.RON Main Menu 



Create a Distance M atrix 

Ihis utility is for the initial entry of a 
matrix into a data file. The user will be asked 
size of the problem. 



the homeport and aids names. 



distance 
for the 
and the 




distances. The size of the problem is the number of aids 
plus one, for homeport. The utility will next ask for names 
of homeport and the aids to be visited. These names will be 
associated with a row and column number of the distance 
matrix. The names are text entries up to a maximum of 15 
characters. 

The utility will next ask for distances. It will 
ask for the distance from a row number to a column number. 
If the user sets up the distance matrix as shown in Figure 
A.1 he will be able to then just read each row across its 
columns while entering the data. These distances must be 
real numbers. The distances should be nautical miles and 
may range from .01 nm to 9999.99 nm. Since these entries 
are real numbers there must be a decimal point in each 
entry. If, for example, the user needs to enter two 
nautical miles, he should enter it as 2. <return>. Be 
careful with these entries; if a decimal point is not 
entered the distance will be off by a factor of 100. For 
quickly and easily entering data the numeric keypad on the 
keyboard may be used and <next> may be used in place of 
<return> to enter the data. If an error is made while 
typing in numbers and <return> or <next> have not yet been 
pressed, then the error may be corrected by using the back- 
space key and retyping the number. If an error is made and 
<return> or <next> have already been pressed, then the user 
should note the row and column number of the error and use 
the cor rec t ion utility to correct the entry. 

When the user finishes typing all entries the 
program will prompt for a filename in which the data will be 
stored. This name may be up to 15 characters long. If for 
some reason the user needs to end the terminal session 
before the entire matrix is entered, data entry may be 
stopped by typing in a negative distance. The user will 
then be asked for a filename for the partial data file. To 



48 



complete the partial data file the user must use the correc- 
tion utility. 

2- Ex t rac t a Ma t rix 

This utility is designed so tnat the user may choose 
a subset of aids frcm a larger matrix. The utility will 
create the proper data file for a subset of aids. 

The user should have an idea of all possible aids 
which may be visited on a general buoytender run (i.e. the 
spring run, the north run, the inside run, etc.) . The user 
should create a distance matrix using the cr ea te utility of 
all aids which could possibly be visited on a general run. 
Usually cnly some portion of these aids will be visited on 
an actual run (for example on a particular spring run only 
90 % of the aids which might be visited on this run need to 
be visited) . In this case the user creates a matrix of all 
aids which might be visited on a run and then each time this 
run is scheduled the user can then use the ext ract utility 
to guickly and easily create a data file for those aids 
which need to be visited. 

Before using the ex tract utility it is suggested the 
user first use the list utility on the larger or initial 
data file. This will provide the user a list of aid names 
in the data file and present them with their corresponding 
row number. 

The ex tr ac t utility will prompt the user for the 
name of the data file containing the general run’s aids. 
The user is then asked for the size of the new problem 
whereby the user enters some m, m less than n in the initial 
data file. The utility then asks the user m times for row 
numbers from the initial matrix. After m row numbers have 
been entered the utility creates the new data file and askes 
the user for a filename for this new data file. 
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3 . C orrec t a Pr e vious ly S^red M atrix 

This utility is to correct entries in a distance 
matrix or tc complete a partially entered data file. The 
utility has a menu which askes if the user wishes tc make 
single entry or sequential entry corrections. The single 
entry mode asks the user for the row number and column 
number of the entry desired to be corrected. The user then 
is asked for the new distance and is returned to the correc- 
tion menu where he may select to make another correction or 
exit the utility. 

The sequential corrections are similar to the entry 
of distances in the c reate utility. The user is asked for 
the starting row number and starting column number where 
corrections are to begin. The utility asks for the 
remaining column entries on the starting row and then will 
proceed to the next row until entries are terminated. The 
utility automatically terminates at the end of the matrix or 
the user may terminate the sequential mode early by entering 
a negative number. If a negative number is entered it will 
not appear in the distance matix. Terminating the sequen- 
tial mode returns the user to the corrections menu for 
another correction or to exit the correction utility. 

Dpon exiting the cor rec tion utility the user is 
asked for a filename for the corrected file. This filename 
may be up to 15 characters long. 

Prepar e a D istance M atrix for Hard cop y 

This utility is for preparing a hardcopy of a 
distance matrix from a stored data file. The data files are 
saved as sequential files, (i.e. , the data files are a 

string of numbers) . This utility will construct a file 
which will resemble a matrix. Due to the width and length 
limitaticns of the printer and paper, the matrix is printed 
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in sections which will have to be cut and pasted together to 
display the whole matrix. A hard copy of the whole matrix 
will provide the user with a means of checking to see if 
data entries are correct. If an entry error is found, the 
matrix also helps the user to locate the row and column 
number of the error. 

Ihe utility asks the user for the name of the data 
file from which a distance matrix is needed. The utility 
creates the hardcopy file and then prompts the user to 
obtain a hardcopy by using the FORMAT command to print out 
file BATFRT.DATA. 

5- Names 

This utility provides a hardcopy of the aid names 
from a data file. This may be useful to the user when using 
the extract utility. The user is asked for the name of the 
data file from which the list shall be made. The utility 
then prepares a file containing the row numbers and their 
associated aid names. The user is prompted to use the 
FORMAT command to obtain a printout of LIST-DATA. 

These matrix utilities are for assisting the user in 
creating and manipulating data files and distance matrices. 
A correct data file is necessary for the proper operation of 
the buoytender route selection program, ROUTE. RUN. 

B. USE CF ROUTE. RON 

The user of ROUTE. RUN is very simple once a data file 
has been created. The user runs the ROUTE. RUN program and 
he is prompted for the name of a data file. The user is 
then asked for the name of an output file for the results. 
The program obtains the size of the problem from the data 
file and will compute an approximate solution time. This 
approximate solution time is based on results from some test 
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problems solved by the program and it assumes the user has 
version 8.01 of the operating system. If another version of 
the operating system is used the solution time may differ. 

After displaying the approximate solution time the 
program reads the rest of the data file and proceeds to do 
the necessary computations. The user, depending on the size 
of the problem, will have just enough time to get a cup of 
coffee before obtaining the solution. The solution will be 
displayed on the screen and will also be placed in the file 
previously named for results, from which a hardcopy may be 
made. The user now has a possible route to take for that 
particular aids to navigation run. 
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APPENDII B 

SAMPLE OOTPOT OF EaOLTENDEH BOOTE SELECTION PfiOGBAM 

The solution for data file BLACKHAW.MAT 
Total tour distance = 455.00 
The tour is as follows 



1 


YB ISLAND 




20 


MOSS LBB MLA 


19 


MONTRY BY 


LB B 


18 


MONTBY HBR 


MB 


17 


PT PINOS LKB2 


16 


PT CYPRESS 


LGB6 


c 


P BLANCAS 


LWB4A 


6 


SAN SIMEON 


LBBl 


7 


VON HELM R 


LGB4 


15 


ESTERO BY 


LWB 


8 


ESTERO GBIOE 


14 


MORRO EY LBBl 


13 


PT BOCHON 


LWB2 


10 


W’ESTDAHL R 


LBBl 


9 


SOUZA R LGB 


12 


LANSING R 


B 


11 


PT SANLUIS 


R B 


4 


SANTA CRUZ 


MB 


3 


AN ISLAND 


LWB8 


2 


MONTAEA LWBIOA 


1 


YB ISLAND 
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APPENMX c 

LISTING OF BDCYTENDEB EOOTE SELECTION PROGRAM 

C MAIN PROGRAM *** PUOYTEKDRP ROUTE SELECTIC'I *** 

IMPLICIT CHAPACTFP*15 (A) 

REAL G, TDIS , TORDIS , DIE ( I 5) 

INTEGER I,J,L,C,P,NUM. ,NOD,SAM(l 5) 

CHARACTEP*15 FN’A''' E , ONA'^ E 
COMMON /NAME/ AMA‘^r(7 0) 

COMMON /WALK/ N,ITOUR{71) 

COMMON /TDST/ TD(7A,7A) 

COMMON /INVT/ IT (14^) 

COMMON /DIST/ 0(70,70) 

COMMON /SAMT/ ISTOR(15,71) 

COMMON /T^.TR/ ITPT0R(1 40) 

CA_LL OBTAIN { FNAM E , ONAM E , M) 

0PEN(3 , F ILE=0''A.ME, STATUS= 'NEW , FORM= ' FORMATTED ' ) 
TORDIS=9999° . 99 
JCALL NODE (N , NUM , SAM) 

CALL N’BR (NUM , SAM , DIS ) 

‘do 7 L=l,5 

DO 8 J=I,N+1 

ITCUR(J)=ISTOR(L,J) 

8 CONTINUE 

CALL MARKDO 
DO 100 C=1,N 

CALL TRPREP(C) 

R = 0 

120 1=1 

IT(1) = I TOUR (1) 

IT (2) =ITOUP (2) 

130 CALL SELCTYd ,C,G, NOD) 

1=1+1 

IF(G. LE. (0 . 0001 ) ) GOTOnO 
CALL ADJTGR (I , R,NOD) 

GOTO130 

110 IF(P. EQ. DGOTOIOO 

CALL '5EVTOP(R) 

G0T012^ 

100 TDIS=0.O 

DO 140 J=1,N 

TDIS=TDIS+D ( I TOUR (J) , I TOUR ( J+1 ) ) 

140 CONTINUE 

IF(TDIS.GF.T0PDIS)G0T07 

T'?'RDIS=TDIS 

DO 5 J=1,M+1 

ISTCP (1 , J) =I TOUR (J) 

6 CONTINUE 

7 CONTINUE 
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VO CO 



r-\ 



20 



29 

S'’ 

31 

32 



23 I = l,\’ + 1 
ITCCR ( I ) =I STOR (1,1) 
CCMTl'.'LT 



CALL TRPPEP(l) 
WPITE(*,22) FMA'-F 
WRIT" (3, 2?) FNAWE 
W R I T R ( * , 3 9 ) T ( 



!IS 



WPITE(3,32) TOFDIS 
WRITE (*,31) 

WRITE (3,31) 

WRITE (* , 32) ( ITCL’P ( I ) 
WRITE (3,32) ( ITCUR ( I ) 
CLOSE (3 ) 

STOP 



, AMAWF ( ITOUR ( I ) ) , 1=1 , W+1 ) 
, ANA^'E ( ITOUP ( I ) ) , 1 = 1 , N'+l ) 



FORMAT ( 1 X ,' Tho solution for data file ’,A) 
FORMAT ( IX, ' Tot il tour distance = ',F8.2) 

FORMAT ( 1 X ,' Th e tour is as follows:') 

FORMAT (IX, 13, 5X, A) 

END 



S UB ROUT I '■ E ORTA I X ( FMA'-’ E , OMAM E , M ) 

IMPLICIT CHARACTER* 15 (A) 

INTEGER N 
REAL TIMF 

CfU" FACTEP* 1 " M ’'lA''' E , ONAW F 
C"''MON /NAWE/ ANAA'E(72) 

CC :M0N /DTST/ D (70,70) 
i-vPITE (* ,5) 

FEAC(*,7) FNAMF 
WRITE (*,’=•) 

P,EAD(*,7) OXA'-'E ■ 

OPEN (4 , FILF = F^’A"E, FORM= ' FORM.’' TTED ' ) 

READ (4 , , FPc = l 4 ) 

TIME=.5"4 + .o2*N + .oqi<l*M**2 + .‘^99R''’994*N**" 

WRITE (*,!"') TIME 
W P I T E ( * , 1 1 ) 

READ (4 ,12,FRF = i 4) ( (D(I,J) ,J = 1,N) ,1 = ] ,N) 

- R "A D ( 4 , 7 , E P R = 1 4 ) ( A N AM E ( I ) , I = 1 , N ) 

CLOSE (4 ) 

P ETL'PN 

14 'W R I T E ( * , 1 5 ) 

CLOSE (4 ) 

STOP 

FOP” AT ( 1 X , ' Un '■ a r v/nat file nare is tn.e fistanc- 
FORMAT ( ■ ) 

FORMAT ( 1 X, ' Ei 1-^ ne.re 
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a t r i y ? ' A ) 



Q 

^ 1'"' 

11 

12 
15 

C 



3 



4 



2 9 



•7 



?0py M (13) 

5’0pv.’\T (1X//IX, ' It v;ill take aapr ox itia te 1 y ',F*5.2,' so 
FORK AT ( 1 X , ' to solve this oroble^. Please '■•’ait.'///) 
FCRKAT (F-. 2) 

F0R;'''-T (IX, ' Fr tor uoon load’) 

END 



SUBROUT P'R. \'CUF (X , r.'UK , SAK.) 

lUTFGEP N , J, K, ITFXP, MUX , SA.'KJ^S) 
CC-'YCM /RAND/ IX,IY,IZ — 

IX=I52o3 

I Y = 2 5 4 R ■ / 

12=32179 

MUr' = 5+M IMT (2 . C*ALOG (FLOAT (^'-4 ) ) ) 

DC 2 J=1,MUM 

I T EM P I MT ( R A ( 1 . '^ ) * F LvOA T ( M ) ) 

IF ( ITEXP. EO.^)GCTQ3 



F=n 

DC 1 K=1,MUM 

IF ( ITEXP. FQ. SAXi(K) ) F = I 

IF ( F. EC. 1) C^T02 " 

SAy ( J) =ITEyp 
CONTIECF 
RETUP'^ 

E n 






y 



/j; 



p >v t 






.Z 7 ' y ^ 

r 



( 



A 









SUEFCUT r:F );r.p (mum , sax , dis ) 

IXTFCER MUX, I,J,K, IMCDF, I TEX P , L , S AX ( 1 - ) 
REAL TFVPU, FIS (1 5) 

COXX^M /MALK/ M,ITCUR(71) 

C 7X.''CV /FIST/ C ( 7 0 , 7 2 ) 

CCX.MON /SAX’T/ ISTCR(I5,71) 

CC'^XCM /TDST/ TD(70,7'"') 

DO IF L=1,MUX 



IMODE=SAX(L) 

DC 20 1=1, M 

DO 2? J=I,M 

TD(I , J) =D (I , J) 
CCMl IMUE 



mcf:^ 



TD (I , IN'CDE) =9999. 99 
C y N T I ' i o ID 

ISTOR (L, 1) = i::ODE 
DIS (L) =E. 0 
DO 21 I = 

I ^ \4 p - 1 



con-5’ ) 



- * I 
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22 



2 ? 
21 



13 

235 



2 2 
210 



C 



A 1 



7pvpp = TD (I3TOF (r. , I) , I FE.MP) 

DC 22 J=2,^; 

IF(TD(irTOR(L, I) ,J) .CF.TEr'PD)G0T0 2 2 
ITDHP=J 

TE''1PD=TD ( I STOP (L , I ) ,J) 

CONTINUE 

ISTOR (L , I +1 ) = I TF^'P 
DIS (L) =D IS (L) +TFNPD 
DO 23 K=1,N 

TO (K , ITF'^P) ='!'nP9 . 

CONTK.'JE 

CONTINUE 

ISTOR (L , N+1 ) = I STOP, (L , 1) 

D IS (L) =D IS (L) +0 ( I STOP (L , M) , ISTOR (L , 1) ) 

CONTINUE 

N0EX=1 

DO 213 I=1,NUM-1 

IF(DIS(I) .F0.DIS(I + 1) )DIS(I+1)=9'^99?.'^9 

IF(DISd) .L'^.DTS(I+1) )G0T0213 

NCEX=3 

TEMPD=D IS ( I -^I) 

DIS (I-i-l) =DIS (I) 

D IS ( I ) =TFXPD 
DO 22'^^ J = 1,N+1 

ITE-'P=I STOP (T +1 , J) 
IST0R(I+1,J)=IST0P(I,J) 

ISTOR (I , J) =ITFO'^P 
CCr'Tn.UE 
CO^'TI^;UE 

IF(NOEX. NT. 1)G0T02"'5 
RETURN 
END 



SUE ROUT I NT TP PREP (C) 

INTEGER G,J,K 
COMMON /WALK/ N',ITCUR(7 1) 
COMMON /TM TP/ I T PT OP ( 1 -i '^ ) 
J=3 

DO 4 K = 1 , 

ITPTCr (K) =ITOfT (K) 

IT PTC? (K-^N) = I TOUR (K) 

IF( ITPTOR fK) . F/.C) ,I=K 
CDNTP'UF 
DC 41 K = l,': 

ITOUR (K) = ITPTOP (J-^"-l ) 
CONTI'UF 
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ITCL’R (N + i ) =I (1 ) 

RFTL'PN' 

KND 



SUBROUTIN'E ''.^RKP {) 

INTFGFR J,K 

COMMON ^LK / N , I TOUP. (71) 

COMMON /TDFT/ TD(7F,77) 

DO 57 J=l,v 

DO 51 :< = 1,N 

TD(J , M) =-l . 7 

51 CON'TINUE 
TD(J, J) =1.0 

50 CONTIM'JE 

DO 52 J = l,';-1 

TD(ITOUP (J) , ITOUR(J + l) )=!."' 

TD (ITOL'P (J +1) , I TOUR (J ) ) =1 . 7 

52 CONTINUE 

TD ( I TOUR (N) , I TOUR ( 1 ) ) =1 . 0 
TD(ITOUP (1 ) , I TOUR (N) ) =1 . 0 
RETURN 
END 

S UB P^UT I N F 5 ELCTY ( I , C , C , NOD) 

INTEGER I,C,'-',F,COE 

IN'TEGER E,ST\'(5) ,RTL5(5) ,'-'D(5) ,NGD 
F FAL G, TEMP, SY (5 ) , SLA ( 5) , STO P , TAP PA Y ( 7 "^ ) 
COMMON /aALK/ N,ITGUP(71) 

C'^r/v.3^: /test/ TD(7 0,'”') 

COMMON /INVT/ I T ( 1 -"l '^ ) 

C OMM.C i: /D I S T / C ( 7 7 , 7 0 ) 

DO ^5 J=1,N 

r; ^RAY (J ) =TD ( ITOUR (2 ) , J) 
o5 CONTI'MJE 

DO SC ;< = 1,5 
F = 0 

T EM P = 0 9 7 ^ on 

DC 61 J=3,N 

IF(D (IT"UR (2) , ITOUP (J) ) . CD . TE'*' F) GOTO'' 1 
IF (TAPRAY ( T TOUR (J) ) . CT . D . '')G7t’0'R1 
TEMP=D (ITOUP (2) , ITOUR (J) ) 

F=J 

61 CONTINUF 

IF ( F. NE. '’) GOT'^-4 
:-'=K-l 
K = 5 
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64 



GCTC6"' 

TAFFAY ( ITQUP (F) ) =l . C 
3Y (K) =TF^^.P 

SLA (K ) =D ( ITOUR ( 1 ) , I TOUR ( 2 ) ) - D ( T TOUR ( F-1 ) , ITFUR. ( 1 ) ) 
STr: (K) = ITCUP (F) 

STLA(K) = I TOUR (F-l) 

S LA ( K ) =S LA ( K ) +D ( I TOUR ( F) , ITCUP ( F-l ) ) -S Y (K) 

UD (K) =F-1 
^ =K 

6C CONTIL'UE 

STOP=n . 0 

IF ( r-’. EO. A) F T0P = 9«'^9 . 90 
TEMP = S LA ( 1 ) 

CCL=1 

DO 63 K=1,M 

I F ( S LA ( K ) . LF . TF'-1 P ) G0T06 3 
TFMP = S la ( K ) 

CGL = K 

63 CCNTIMUE 

G = 3LA ( COL) +STCP 
^:CD = ND (COL) 

IT (2 * I+l ) =STF (CCD 
IT (2* I +2) =STLA (CRL) 

RETURN’ 

END 

r 

SUBPCUT IMF ADJTOR ( I , R, NCD) 

IMTECEP I,J,R/MCO 
CCMMCM /WALK/ M, ITOUR (7 1) 

COiD^M /TOST/ TD(7''’',72) 

CCM'^S'N /IMVT/ IT (140) 

COMMON /TMTP/ ITPTOR(14r') 

OC '=1 J = 2,NCD 

ITPTCR (MOD + 2-J) = I TOUR (J) 

31 CONTINUE 

DO 82 J = 2,iIGD 

ITOUR (J)=ITPTOP(J) 

3 2 COST IN UK 

TO (IT (2* I -2) ,I'7(2’"I-1) )=1. '' 

TD(IT(2*I-1) ,IT(2*I-2) ) = 1D 
TD(IT(2*I-1),IT(2*I))=1.R 
TC(IT(2*I) ,IT(2*I-l) )=1.0 
TE(IT(2*I) ,IT(1))=1."’ 

R = 2 

PF.TUPN 

END 
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SUoFOUT I VE PEVTfF(P) 

I\'TFGER P 

COMMON /.%'^LK/ G,ITOUP(~l) 

co'/yor. /T'”rn/ 

DO J = I,\'+1 

ITPTOF (\’ + 2-J) =I'"OUF (J) 

9? CCKTIKLIE 

DO 91 J=] , N+1 

ITCO? (J) =IT?TCP (J) 

91 CO.NTI'OJE 

R = 1 

RETL’P^' 

END 

EENCTIOM RA-'OC^CA) 

REAL TFO.F 

CO^’,MO^: /RP'jn/ Iv,iY,IZ 
I X=171*’10C ( I X, 177) - 2* ( I X/1^7) 

I Y=172*XCD ( I Y, 1 7E) -2* ( I Y/1 76) 

IZ = 170*YOC(IE ,17?) - 2* (I7,/17 = ) 

IF(IX.LT.O) IX=IX+3?269 
IF( lY. LT.^) IY=I Y+3^3?7 
IF(IZ.LT.'")I7 = IZ + 37'’-23 

T EX P = F LCAT ( I X ) /;^ 9 2 6 9 . +F L<^-AT ( I Y ) /3 9 3 +F L-9 T ( I 7 ) /"’ '' 3 2 ’ . 

R A XDCM = A CD ( T EY r , 1 . 9 ) 

R FTUR Y 
EL'O 
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APPENDIX D 

LISTING OF MATRIX UTILITIES PROGRAM 






12 



1 

2 



5 

5 

7 

2 

9 

19 

11 



*** DISTAMCE MATRIX UTILITIES *** 
IMPLICIT Ch'A RACTER* 1 5 (A) 

IMTECFR I A MS 

COt-M^.ON /DIST/ D (l''*'^, 10 9) 

COMMON /NAME/ AMAMF(1'''’) 

WRITE (* , 1) 

WRITE (* , 2) 

WRITE (* , 3) 

WRITE (* ,4) 

WRITE ( * , 5) 

WRITE(*,6) 

WRITE (* ,7) 

WRITE (* ,8) 

WRITE (* ,9) 

WRITE (*, 19) 

READ (*,11) IANS 

IF ( ( I ANS . LT. 1) .CR. ( I -MS . GT. 6) ) GOTO 12 
IF( IA^:S. FQ. 1)CALL "ATRIXO 
IF( IANS. EC. 2)CALL EXTRCT() 

IF ( IANS . EQ. 3)CALL CORECT() 

IF(IAMS. E?.4)CALL MATPRTO ' 

IE( IANS. FQ. 5)CALL LIST() 

IF( lANS.NE. 6)CCT012 
STOP 



F 9RMAT ( 1 X , ' Th 1 s is a set of utilities to be 
F ■ RM.'-.T ( 1 X , ' the orootarr ROUTE. RUN. ,^hsy ire 
EOR.MAT ( 1 X , ' c rea te an<4 nanipulate a 'iistance 
1 - create a Iistance -natrix') 
- extract a sr^all Tatrix frc"', 



FORMAT (//IX, 
FORMAT (IX, ' 2 
FORMAT ( IX , ' 3 
FORMAT (IX, ' 4 
FOPM'-T (IX, ' S 
FORMAT (IX, ' 5 
FORMAT (//IX, 
FORMAT (I 1) 
END 



used with') 
to help ' ) 
matrix.' ) 



^ laraer aatrix') 

- correct a nroviously stored aatrix')^ 

- create a hardcopy of a distance nntrix') 

- list naiTies with their row nuaiber in matrix') 

- exit this proorar.') 

''.■'hat is your selection? '\) 



SUBROUTINE MATE IX () 

IVPLI-2IT CHA RACTFR* 1 R (A) 
REAL DIST 
INTEGER N, I , J, L 
C0^''^’0^: /DIET/ D(19'^',1^9) 
COMMON.' /NAME/ A ’':A F ( 1 ''' 9 1 
W R I T E ( * , 3 1 1 

READ (*,32) N 
WRITE (*,33) 
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DO 34 I=1,K 

WRITE (*,3 5) I 



READ ( 


*,36) AMAME(I) 


COr.’TINUE 




WRITE ( * , 


37) N,N 


WRITE (* , 


38) 


DO 39 1= 


1,N-1 



D ( I , I) =9<;>99 . 99 
DO 40 J=I+I,N 

WRITE (*,41) I,J 
READ (*,42) DIST 
D(I,J)=DTST 
DO 10^ L = l,100'’ 



100 CCN'TIt.'UF 

D (J , I) =D 1ST 
IF(DIST. GF. 0)GCTQ40 
J 

I=N-1 

41 CONTINUE 

3 9 CON’TINUE 

D(N,N)=99?9.99 
CALL STORE (N) 

RETURN 

31 FORNAT (//IX, ' How larqe is the ocoble"’? '\) 

32 FORMAT (BN, 13) 

33 FORM.AT (//I X , ' I n put names for homeport and aids') 

35 FCF'’AT(1X, 'NuTber' , 13, ' Name? ' \) 

35 FORMAT (A) 

3 7 F0R^:AT(1X, 13, ' BY ',13,' Matcix') 

39 F0PMAT(' Input 'distances from point I to point J') 

41 FOPMAT(' Distance from pt ',13,' to nt ',13,' = '\) 

4 2 rOR:’AT(BN,F7. 2) 

EN'D 



SUBROUTINE EXTRCT() 

IMPLICIT CHARZ^.CTEP* 1 5 (A) 

I NTEGEP N , M , I , J , I ROW (100) 
PEAL SD(100,iap) 
CHARACTFP*15 ASMMM(1'"'^) 
COMMON /DIET/ D ( 1 ^ 0 , 1 '^ ■'^ ) 
Cr*-'MO>: /NAME/ '•NAME(1^'’) 
WRITE (* , r~') 

WRITE (* ,11) 

WPITE(*,12) 

'WRITE (*, 13) 

CALL OBTAIN (N) 
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21 

22 

2 0 

24 

23 

10 

11 

1 2 

1 

14 

15 

C 

16 



i'.-RITE { * , ] 4) 

I' r i D ( * , 1 5 ) 

DO 20 1=1, M 

WRITE (*,21) r 

FORMAT ( 1 X , I 3 , ' . Rov; nu^b^i: Iron Tiain -".atrix: ’\) 

RFAD(* , 15) irOW( I ) 

DO 22 J=1,M 

SD ( I , J) =D ( IROW( I ) , J) 

COMTIM'JF 

ASM MM ( I ) =A XAM F ( I ROW ( I ) ) 

COXTIKUE 
DO 2^ J^l,M 

DO 24 I=1,M 

D ( I , J) = S D ( I , I'^OW ( J ) ) 

COMT I ''.’UE 
D (J , J) =90^'^. 09 
AMXME (J) = ASMX:-' ( J) 

CC.MTINOr 
CALL STCRFC'D 
RETURN’ 

FORMAT ( //I X , ' To create a sr.aller natrix frO:T your nain') 
FCRM.9T (lx, ' ] istance "atrix, enter tF ' row nuToers of) 
Fop\< a T ( 1 V , ' aa i n "^atrix ails which yo.i Mesire in the ') 

F C ? • ’ A T ( 1 X , ' s~ a 1 1 0 r r a t r i x ' ) 

FORMAT (//) X, ' ;.'h at is the size of the smaller matrix? ' \) 
FORM AT ( DU , I 3 ) 

END 

SUFPCUTrX’E CDRECTO 

IMPLICIT CHARACTFR* 1 5 (A) 

INTFGF? '1, I, J,K, I-'C:S, IP, IC 
F. E L D I S T 

C C >' M C N /D I S T/ D ( 1 , 1 ^ ^ ) 

CCMMOi: /f’.AME/ A N’A^’ F ( 1 '^ '’' ) 

CALL OBT.’\IN(t:) 

WRITE (*,14) 

WRITE ( * , I'^) 

WRITE (*, 17) 

WRITE ( * , 1 ^ ) 

READ (*,19) I'.N’S 
IF( IAN’S . F9. 3)CCTC35 
IF ( r ^I’S . FC. 2)rOT''C5 
I F ( IANS . NF . 1 ) GOTO] 6 
WRITE (* ,21) 

P EA D ( * , 1 2 ) IP 
w P I T e: ( * , 2 : "I 
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READ (*,12) IC 

WRTTr. (*,2 1) IP, IC,D(IP, IC) 
vJRITE (* , 24 ) 

RF,AD(*,13) DIET 
D(IR, IC) =DISC 
DO IP''' :< = i,iopr: 

IPO CCIITINUE 

D ( IC , IR) =D 1ST 
G0T016 

25 WRITE(*,2P) 

WRITE (* ,27) 

READ(*,L2) IP 
WRITE (* ,23) 

PEAD(*,L2) IC 
DC 3 0 I = IR,:;-L 
DC 31 J=IC,M 

WPITF(*,2^) I,J,D(I,J) 
v:rite (* ,24) 

READ(*,13) dist 
IF(DIST.GE. P. '"I) GOTO 3 2 
I =N 
J=\: 

GCTC31 



32 


D (I , J) =DIST 








DO IIP K = 1 , 1 P^^P 






IIP 


CCN’T INL’C 








D (J , I) =DIST 






31 


cowT r:uE 








IC=I +2 








D ( I , I) =<3^99 . 9 ^ 








CCCTIN’LT 








GOTO 15 






35 


CALL STORE (!') 








RETURN 






12 


FORMAT (3N, 13) 






1 3 


FORMAT (F7. 2) 






1 A 


FORMAT (//I X,' Th i s utility is foe oocre:' 


tine or corolatin^' 


) 


15 


FORMAT (IX, '3 iroviously storoF Psta fil 


e . • ) 




1 7 


FP'RM AT ( 1 V , ' 1 - sip.ole -^ntrios 2 - 


Seouenti al entries' 


) 


1? 


F0RMAT(1X,'3 - Exit corrections 


5-c lection: ' \ ) 




19 


FORMAT (11) 






21 


FCTMAT (lx, ' Row nir-^ber to correct ' \) 






22 


F or;-; AT ( 1 X ,' Go 1 uci n na^nber to correct 'A) 






2 3 


FCPMAT(1 X, ' Pres_-nt Oistance froT,',I3,' 


to ’ , 13, ' = ■ ,-7. 2) 




24 


FOr;‘.AT(lX, ’ -'cw list. a nee = ' \) 






2 6 


FC^'-'AT (1 V , ' To en-'' .antry of Pata, f/ne -a 


n ?o a t i e n i s t ~ n c r ' 


) 
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27 

2P 

7 

o 

-/ 

IG 

3 

3 

4 

5 

5 

11 

1 2 

1 3 



F ORy.AT ( 1 y , ' S t a r t i ng rov/ nu~' ber ' \) 

FORFAT ( ly, ' Star ting coluTn n j.rbcr '\) 

EKD 

SUBPOUTI'JG yATPRTO 

IMPLICIT CHARACTFR*1F (A) 

I\'TFGFR I, J, y, IS, IF, JS , JE 
COMyON' /FIST/ D(1^'’,1GG) 
coyyoM /MAr.F/ Ai’AyE(i"iG) 

WPITE (* , 3) 

WRITE (* ,4) 

WRITE (* , 5) 

WRITE (* ,6) 

CALL CETAI,‘‘'(M) 

OPEN (3 , F ILE= 'y.ATPPT. DATA ' , STATUS = ' \’ EW ' , FORy= ' FCRy.ATTEE ' ) 

I S = 1 
I F=M 

IF (U , GT, 5'^) I ^^ = 5^ 

JS = 1 
JE = 9 

DO ? I=IS,I“ 

IF(W.LT.1'^) JF=K 
DO IF J=JS,JF 

W R I T E ( 3 , 1 1 ) D { I , J ) 

C ONT I 'JF 
WRITE (3,12) 

CO^iTIEUE 
WRITE (2,1') 

JS=JC+1 

JF=JF+9 

IF(JE.CT.C:) JF.=y 
IF(JS.LF.'') GCTO'^ 

I S = I E + 1 
I E = I E + 5'’ 

IF(IE,GT.W) IF=v 
IF(IS.LE.N) GOTO" 

CLOSE (3 ) 

WRITE {* , 14) 

RETURW 

FORMAT (//I X, ' T'-'. is utility orapaces a f i 1 containing tho ’) 
FORMAT (IX, •'=' istance T.atrix for the orint^^ r. You --a'/ 'v^.va') 
FORMAT ( lx ,' to cut 5. oa s 1 3 so-^e sections to net t"^.- -atrix') 
FORMAT ( 1 X. , • i n n by n ' fo re . ' ) 

FORMAT (F'^ . 2\) 

FOF'MAT(' .') 

FORM .AT ( ' ' ) 
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14 



FORMAT ( ' Use FORMAT coa, nand to ootain orint of 'OATRPT.OA 

E MD 

subroutp:e r.isTO 

IMPLICIT CHAPACTEr*15 (A) 

IN]TECER U,l 

CCMMCC /CIST/ D (1 1 PC) 

C OMMON /came/ a ^’A.M F ( 1 S ^ ) 

WRITF (* , 1^) 

WRITE (*, 11) 

WRITE (* ,12) 

WR ITF (*,13) 

CALL ORTAIM(N) 

OPEN (?, FILF = ' L 1ST. nATA ', STATUS = ' , FOPy= ’FCRMATTF" ' ) 

DO 14 1 = 1, :■ 





WRITE (*,15) I, ANA ‘I 


F(I ) 








WRITE (3, 15) I,A.NA*^ 


E(I) 






14 


EONTIi;UF 










CLOSE (3 ) 










RETURN 








10 


F0RMAT(//1X, ' Tnis routine 


will 


list ■’iO row number an' 


11 


FCo^'AT ( IX, ' a i^ nare. 


The 


list 


v.'ill scroll by the ' ) 


12 


FORMAT ( 1 X ,' sc reen ana 


b e 


saveT 


to T i s !< u n T a r file ' ) 


13 


FORMAT (IX, ' LIST. DATA 


U5G 


FORMAT 


coT.Tana for hardcopy' ) 


1 5 


FORMAT (IX, I 3, 5X, A) 








O'* 


END 










SUP 'ROUTINE ■'PT'^IN(l') 









IMPLICIT CHAPACTFR*15 (A) 

INTEGER .\',I,J 
C !•' A F ' CT E F* 1 ^ " v,- M E 
COMC'CN /DIST/ D (1 1'’^) 

COMMON /'-iME/ ANAMF(10"’) 

WRIT- (* , 1^) 

READ (*,11) FMA^F 

OPEN ( 4 , F I LF=F "'A.''* F , FORM= ' F OR^ ' TTFF ' ) 

READ(4 , 12, rpp = i.'i) 

read (4 , 13, E?F = 14) ( (0 (I , J) , J = l, N) , 1 = 1, A') 

READ(4 , 1 1, EF^ = 14) (A‘:A''E(T ) , 1 = 1 , N) 

CLOSE (4 ) 

RETURN 

14 WRITE (*,15) 

CLOSE (4 ) 

STOP 

1'’ FCR'-MAT ( //I X , ’ Un-i ?r '.•/nat file is tWe ^istanca ratrix? ’\) 

11 FORM AT ( A) 
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12 rGDMAT(I3) 

'> TT T7 ■T’ ^ 2 ) 

15 5’G?'-' AT ( 1 X , ' Cr c o c lucinn loo"'.') 

EKD 



SUTPOUTr’f^ STOR'~{X) 

IMPLICIT CHAFACTFP* 1 5 (A) 

r.’TFGCF ",l,J 

CHA?ACTFR*15 
CC''m' 0N /DIST/'D'dPP, 

CO’^IMCX AXA^’Fd'd 

WRITE (*, 10) 

REAP (*,11) FN’A’^F 
WFITE(*,12) FXAr-F; 

0PF'’(3 , FILE=FVAMF, STAT'JF= 'MEW ,FOPM= 'FOP'^ATTCO' ) 
WriTE (3 , 1^, FPF = 16) M 

WRITE (? , 1^ , EP.P = 1 F) ( (D (I , J) , J=1 , X) , I = I , X) 

WRITE (3 , 11, FFF=1F) (AMAME(I ) , 1 = 1, M ) 

CLOSE (3 ) 

WPITE (* , ] 5) 

G GTC 1 P 



16 


i.'RITE { * , 17) 




1 r 


CL<^SE (3 ) 






R ETURN’ 




1'' 


FOPM^T (//IX, ' Fi le mrs for ’latrix to 


11 


FCPFAT ( A) 




12 


f^OR:;AT(' Tha file will be sov 


un<ier 


1 3 


p C ^ V 7 1 1 3 ) 




lA 


FO^'TAT (F7. 2) 




1 5 


^CF” IT (1 X, ' Fi le saved without 


•? r r o r ' ) 


1 7 


for;; AT ( 1 ' Er ro r .vhile saving 


f ila’ ) 



EGG 
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